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CHAPTER  I 


INTRODUCTION 

This  report  describes  the  results  of  a  fifteen-month  study,  sup¬ 
ported  by  DCA  Contract  100-82-C-0026,  of  the  real-time  implementation  of 
an  algorithm  combining  time-domain  harmonic  scaling  (TDHS)  and  Adaptive 
Residual  Coding  (ARC)  at  a  transmission  bit  rate  of  16  kb/s.  The  modi¬ 
fications  of  this  encoding  algorithm  as  originally  presented  by  Melsa 
and  Pande  (1981)  to  allow  real-time  implementation  are  described  in 
detail.  A  non  real-time  FORTRAN  simulation  using  a  sixteen-bit  word 
length  was  developed  and  tested  to  establish  feasibility.  The  hardware 
implementation  of  a  full-duplex,  real-time  system  has  demonstrated  that 
this  algorithm  is  capable  of  producing  toll  quality  speech  digitization. 

This  report  has  been  divided  into  two  volumes.  The  first  volume 
discusses  the  algorithm  modifications  and  FORTRAN  simulation.  The 
details  of  the  hardware  implementation,  schematics  for  the  system  and 
operating  instructions  are  included  in  Volume  2  of  this  final  report. 

The  minimum  requirements  for  the  TDHS- ARC  speech  digitization  tech¬ 
nique  are  that  it  should  provide  toll  quality  speech  under  optimum  con¬ 
ditions  and  high  quality  robust  speech  in  the  presence  of  60  dB  acoustic 
background  noise  and  in  the  presence  of  random  channel  error  rates  up  to 
0.001.  The  algorithm  should  be  designed  for  a  data  rate  of  16  kb/s. 
[Statement-of-work] 

The  trade-off  between  transmission  quality  and  hardware  complexity 
was  studied  with  the  overall  objective  of  finding  algorithm  refinements 
that  simplify  implementation  without  degrading  the  previous  performance. 


It  is  important  to  note  that  minor  sacrifices  in  system  performance  may 
be  made  to  achieve  significant  hardware  simplification.  This  volume 
presents  the  software  simulation  for  a  time  domain  harmonic  scaling  - 
adaptive  residual  coder  system.  It  includes  the  finite  word  length 
simulations  as  the  final  form. 

A  necessary  condition  for  understanding  the  direction  and  signifi¬ 
cance  of  this  work  is  a  thorough  familiarity  with  the  original  form  of 
the  TDHS-ARC  9.6  kb/s  system.  Therefore,  a  brief  survey  of  this  system 
comprises  the  second  chapter.  Under  the  heading  of  "Final  System",  the 
third  chapter  outlines  the  final  algorithm.  This  chapter  will  highlight 
the  major  modifications  without  providing  the  background  associated  with 
each  variation.  The  fourth  chapter  presents  the  algorithm  changes  which 
were  primarily  due  to  the  system  options  called  for  in  the  statement  of 
work.  The  changes  dictated  by  hardware  considerations  are  discussed  in 
the  fifth  chapter.  The  sixth  chapter  reports  the  effects  of  replacing 
the  error-free  channel  approximation  with  a  noisy  channel.  The  seventh 
chapter  is  concerned  with  an  explanation  of  the  program  options  and  out¬ 
put  files,  as  well  as  their  use.  The  final  chapter  forms  a  summary  of 
the  work  and  the  final  system. 
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CHAPTER  II 


ORIGINAL  SYSTEM 


2.1  Introduction 

Throughout  this  report  references  will  be  made  to  the  “original 
system”.  This  is  a  reference  to  the  9.6  kb/s  TDHS-ARC  system  presented 
by  Pande  [Melsa  and  Pande,  1981].  A  block  diagram  which  reflects  the 
major  points  of  Interest  of  the  system  on  this  date  is  illustrated  in 
Fig.  2.1.1,  entitled  “Original  TDHS-ARC  System”.  This  is  essentially 
the  beginning  point  of  the  research  herein  reported.  In  order  to  eval¬ 
uate  the  significance  of  this  work,  it  is  necessary  to  understand  this 
original  system.  This  chapter  also  highlights  the  characterstlcs  of  the 
original  system  which  were  modified  to  obtain  the  final  system;  these 
modifications  are  the  subjects  of  subsequent  chapters. 

The  basic  structure  of  the  original  system  is  illustrated  in  Fig. 
2.1.2.  This  structure  was  first  proposed  by  Malah  in  1980  when  he  pre¬ 
sented  the  results  of  a  TDHS-CVSD  system  and  suggested  the  usefulness  of 
combining  TDHS  with  other  waveform  coders  [Malah,  1980].  The  structure 
is  a  time  domain  harmonic  compression  unit  concatenated  with  an  adaptive 
residual  encoder  in  the  transmitter,  and  an  adaptive  residual  decoder 
concatenated  with  a  time  domain  harmonic  expansion  unit  in  the  receiver. 
The  receiver's  components  perform  the  inverse  function  of  the  transmit¬ 
ter's  sections  (logically  resulting  in  receiver  expressions  such  as  in¬ 
verse  ARC).  The  <}dality  is  limited  by  noise  sources,  such  as,  quantiza¬ 
tion  in  the  ARC  and  compression  versus  expansion  differences  in  TDHS. 
Again,  it  is  worthwhile  to  note  that  this  structure  is  the  heart  of  the 
system  and  was  maintained  throughout  the  investigations. 
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Figure  2.1.1  Original  TD1IS-ARC  System 


The  many  notable  characteristics  of  the  original  system  are  Identi¬ 
fiable  as  pertaining  to  the  transmitter,  channel,  or  receiver.  They 
have  been  individually  classified  as  such  and  are  presented  in  the  fol¬ 
lowing  sections.  Many  of  the  points  of  interest  of  the  receiver  are 
similar  to  those  of  the  transmitter.  They  are,  therefore,  noted  as 
being  congruent  and  only  briefly  mentioned. 

2.2  Transmitter 

The  Input  to  the  transmitter,  in  Fig.  2.1.1  is  denoted  s(k),  an 
analog  speech  signal  sampled  at  6400  Hz.  This  sequence  forms  the  input 
to  both  the  compression  block  and  the  block  pitch  extraction  unit  which 
jointly  form  the  TDHC  algorithm  in  the  transmitter.  A  compression  ratio 
of  2:1  is  used  and  although  many  acceptable*  windows  may  be  used,  a  tri¬ 
angular  window  was  suggested  due  to  computational  simplicity  [Melsa  and 
Pande,  1981].  A  higher  compression  ratio  can  be  effective.  However, 
the  choice  of  a  higher  scaling  factor  entails  the  need  for  a  more  com¬ 
plex  window  function  for  equivalent  performance  [Malah,  1981].  The  new 
characterization  of  block  pitch  extraction  was  precipitated  by  a  major 
topic  of  Investigation  which  Is  described  in  Chapter  V.  A  pitch  period 
la  determined  for  a  block  of  speech  by  locating  the  maximum  of  a  corre¬ 
lation  computation.  The  block  length  and  range  of  the  computation  are 
100  and  20  to  100  respectively.  The  output  of  the  TDHE  unit  is  a  com¬ 
pressed  3200  samples /second  speech  signal. 

*  An  acceptable  window  is  defined  as  any  window  h(i :Np)  such  that 

h(l:Np)  -  1 
h(NpjHp)  -  0 

This  is  required  by  continuity  considerations. 


The  compressed  signal  constitutes  the  input  to  the  transmitter  ARC. 
The  encoder  embodies  a  fourth  order  adaptive  predictor  and  a  21  level 
buffer  controlled  quantiser.  The  adaptation  of  the  sequentially  adaptive 
linear  predictor  coefficients  is  based  upon  the  steepest  descent  of  the 
energy  in  the  error  between  the  predicted  and  actual  values  [Melsa  and 
Pande,  1981].  The  buffer  control  signal  shifts  the  quantizer  thresholds 
to  reduce  the  possibility  of  buffer  underflow  or  overflow.  This  is 
achieved  by  forcing  or  restricting  the  occurrence  of  a  run  length  code¬ 
word.  The  sequence  of  quantization  levels  is  the  only  information  passed 
on  to  the  remainder  of  the  transmitter,  the  channel,  and  eventually,  the 
receiver.  The  receiver  must,  therefore,  require  no  more  than  this  infor¬ 
mation. 

The  coding  block  accepts  the  quantization  levels  and  maps  them  to  a 
variable  length  codeword  sequence  that  adds  the  error  protection.  The 
variable  length  code  precipitates  the  need  for  a  buffered  output  to  the 
channel.  The  buffer  serves  a  dual  role.  It  allows  the  transmitter  to 
emit  bits  at  a  constant  rate  and,  as  illustrated  by  the  buffer  control 
signal,  it  presents  the  adaptation  information  to  the  variable  quantizer. 
The  buffer  outputs  bits  to  the  channel  at  a  transmission  rate  of  9600 
bits  per  second. 

2.3  Channel 

The  original  system  Incorporated  a  noisy  channel  simulation  with  an 
Interactive  bit  error  rate.  The  channel  model  uses  a  pseudo-random  num¬ 
ber  generator  to  determine  the  location  of  a  bit  error  in  the  channel 
sequence .  Given  the  decision  of  the  pseudo-random  number  generator,  the 
algorithm  changes  the  bit  stream  accordingly. 


The  channel  is  constrained  to  a  9600  bits/second  rate.  Although 
this  is  a  channel  constraint,  its  influence  is  accommodated  in  the  re¬ 
ceiver  and  transmitter  algorithms  by  restricting  the  average  bits  allowed 
per  quantization  level.  In  this  way,  it  affects  the  coding  and  decoding 
blocks  and  the  fineness  of  quantization  in  both  ARC  algorithms. 

Therefore,  the  importance  of  the  simulated  channel  lies  in  its 
ability  to  induce  errors  in  the  bit  stream  between  the  transmitter  and 
receiver.  The  transmission  rate  of  9600  bits/second  plays  a  role  in  the 
transmitter  and  receiver  algorithms  and  not  in  the  channel  algorithm. 

2.4  Receiver 

A  bit  stream  at  the  transmission  rate  of  9600  bit/second  is  received 
from  the  channel  by  the  decoder.  The  decoding  block  maps  the  bit  stream 
back  into  a  sequence  of  quantization  levels  which  may  or  may  not  deviate 
from  the  original  quantizer  level  sequence.  The  efficiency  of  mapping 
q(k)  to  q(k)  is  dependent  only  upon  the  specified  error  correcting  capa¬ 
bility,  bit  error  rate,  and  code.  Again,  the  requirement  of  a  sample 
buffer  is  generated  by  the  utilization  of  a  variable  length  code.  The 
sample  buffer  contains  the  quantization  levels  contaminated  by  noise 
induced  in  the  channel. 

These  levels  are  accessed  by  the  inverse  ARC  block.  This  ARC  is 
also  a  fourth  order  adaptive  linear  predictor  combined  with  a  21  level 
quantizer.  Note  that  there  is  no  need  for  a  buffer  control  signal  in  the 
receiver  and  that  the  receiver  need  not  know  the  threshold  values  used  by 
the  transmitter.  Other  than  these  major  disparities,  the  inverse  ARC 
operates  as  the  exact  inverse  to  the  encoder  in  the  transmitter. 
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Th«  time  domain  harmonic  expansion  element  in  the  receiver  acts,  as 
stated  earlier,  much  the  same  as  the  compression  blocks.  A  pitch  extrac 
tlon  routine  Identical  to  that  in  the  transmitter,  is  exploited  in  the 
receiver.  The  major  relevant  differences  are  the  expansion  ratio,  which 
is,  of  course,  1:2,  the  exact  Inverse  of  the  compression  ratio,  and  the 
window.  Here,  the  choice  of  windows  is  forced  due  to  the  choice  of  a 
triangular  window  in  the  transmitter.  Once  again  employed  is  the  notion 
of  an  acceptable  window.  That  is,  one  that  maintains  the  continuity  at 
its  fringes,  as  well  as  those  of  the  succeeding  window.  The  impelling 
concept  of  this  algorithm,  as  a  whole,  is  the  production  of  a  high 
quality  sequence  of  estimates,  s(k),  to  the  original  speech  sequence, 
s(k).  This  estimate  forms  the  output  of  the  receiver  and  hence  the 
algorithm.  It  is  considered  a  6400  Hz  sampled  speech  signal  estimate. 


CHAPTER  III 


FINAL  SYSTEM 

3.1  Introduction 

A  block  diagram  of  tha  final  systea  la  Illustrated  In  Fig.  3.1.1. 
This  simulation  Incorporates  the  specifications  cited  In  Chapter  I.  It 
Is  laaedlately  apparent  that  the  final  algorithm  maintained  the  basic 
structure  discussed  in  the  opening  remarks  of  Chapter  II  and  Illustrated 
In  Fig.  2.1.2.  The  transmitter  preserves  the  TDHC-ARC  format  aud  the 
receiver  upholds  the  inverse  ARC-TDHE  structure.  Given  that  the  major 
concept  of  a  TDHS  -ARC  system  is  maintained,  the  points  of  interest  of 
the  final  algorithm  are  the  modification  of  the  Individual  components. 
Although  detailed  discussions  of  these  changes  are  deferred  for  future 
chapters,  they  are  highlighted  in  the  remainder  of  this  chapter. 

3.2  Transmitter 

The  Input  to  the  transmitter  is  composed  of  an  analog  speech  signal 
sampled  at  8000  Hz.  This  signal  forms  the  Input  to  both  the  compression 
block  and  the  pitch  extraction  block.  The  compression  ratio  of  2:1,  is 
preserved  In  this  final  algorithm.  Also  preserved  Is  the  utilization  of 
a  triangular  window,  the  length  of  which  Is  twice  the  detected  pitch 
period  length. 

The  pitch  extraction  block  underwent  crucial  transformations.  The 
Input  sampled  speech  is  first  low-pass  filtered  by  a  9th  order  linear 
phase  FIR  digital  filter  before  pitch  detection  is  performed.  Following 
this  operation,  a  serial  pitch  extraction  routine  is  executed.  The 
investigation  of  this  routine  constituted  the  majority  of  the  effort 
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Figure  3.1.1  Final  TDHS-ARC  System 


allotted  to  the  modification  of  the  TDHS  algorithm.  It  is  based  on 
hardware  considerations  presented  in  Chapter  V.  Pseudo-random  pitch 
picking  and  a  new  processing  of  the  correlation  function's  data  are 
exercised  by  the  changes.  The  hardware  simplicity  associated  with  the 
three  value  center  clipped  autocorrelation  computation  was  retained. 

The  adaptive  residual  encoder  found  in  the  transmitter  incorporat¬ 
ed  a  fifteen  level  fixed  quantizer  and  a  fourth  order  linear  predictor. 
The  adaptation  procedure  is  modified  to  eliminate  the  division  opera¬ 
tion  called  for  in  the  original  system.  This  is  accomplished  by  adopt¬ 
ing  a  fixed  gain  adaptation  procedure.  Another  division  operation  is 
eliminated  by  using  the  quantizer  normalization  factor  as  a  multiplica¬ 
tive  adaptation  to  the  thresholds  and  scaling  factors  as  opposed  to  a 
dividing  adaptation  to  the  error  signal.  This  results  in  only  a  mild 
truncation  of  the  original  system's  performance.  The  division  opera¬ 
tions  are  eliminated  to  make  the  algorithm  hardware  compatible  with  the 
NEC7720  signal  processing  chip  which  does  not  support  a  division  opera¬ 
tion.  The  use  of  a  fixed  length  code,  a  sampling  rate  of  8000  Hz  and 
the  restriction  to  a  16  kilobit/second  transmission  rate  dictated  that 
the  maximum  number  of  quantization  levels  be  15. 

Due  to  algorithm  changes  presented  in  Chapter  IV,  the  quantization 
levels  became  the  information  passed  from  the  transmitter  to  the  receiv¬ 
er.  The  15  quantization  levels  require  4  transmission  bits/sample  re¬ 
sulting  in  an  effective  transmission  rate  of  4000  quantization  level 
samples/second.  The  quantization  levels,  therefore  formed  the  output  of 
the  transmitter. 
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3.3  Channel 


A  temporary  approximation  to  a  noiseless  channel  is  instrumental  in 
the  reduction  of  the  transmitter  and  receiver  algorithms  to  their  form 
in  Fig.  3.1.1.  The  performance  of  the  system  in  a  noisy  channel  is  re¬ 
ported  in  Chapter  VI.  A  crucial  consequence  of  the  noiseless  channel 
approximation  is  that  the  output  of  the  transmitter  is  exactly  the  input 
to  the  receiver.  In  essence,  the  channel  is  not  present,  save  the 
transmission  rate  characteristic. 

The  channel  transmission  rate  of  16  kilobits /second  is  another  of 
the  channel's  features  that  dictated  significant  alterations  to  the 
original  system.  Although  these  variations  did  not  predominate  in  the 
channel  model  itself,  they  are  critical  to  the  final  algorithm.  The 
transmission  rate  restricted  the  quantization  of  the  ARC  and  thus, 
upper-bounded  the  signal-to-noise  ratio  of  the  ARC  system.  The  effects 
of  transmission  rate,  then,  are  sustained  in  the  transmitter  and  receiv¬ 
er  and  not  in  the  channel. 

3.4  Receiver 

The  channel  output,  which  exactly  coincides  with  its  input  and  the 
transmitter  output,  becomes  the  input  to  the  receiver.  It  is  considered 
a  16  kilobit/second  signal  representing  4  bit  codewords  of  the  quantiza¬ 
tion  levels.  Since  the  quantizer  levels  form  the  transferred  informa¬ 
tion,  no  decoding  is  necessary  and  inverse  ARC  is  applied  directly  to 
the  input.  Inverse  ARC  in  this  system  is  equivalent  to  the  feedback 
loop  of  the  ARC  in  the  transmitter.  It  contains  a  fourth  order  linear 
predictor  as  in  the  original  system.  The  quantizer  has  15  levels  and  is 


fixed  In  the  final  algorithm  as  opposed  to  the  21  adaptive  levels  used 
in  the  original.  The  attractive  characteristic  of  the  inverse  ARC  (not 
requiring  the  threshold  values )  is  maintained  in  the  final  algorithm. 

The  compressed  signal  y(k)  is  a  4000  Hz  sequence.  It  is  the  esti¬ 
mate  of  the  compressed  speech  signal  in  the  transmitter,  y(k),  and  is 
the  driving  sequence  of  the  expansion  and  pitch  extraction  units  of  the 
receiver.  The  expansion  block  of  the  receiver  behaves  exactly  as  that 
of  the  original  system.  This  block  perpetuates  the  1:2  expansion  ratio 
and  the  triangular  window  suggested  by  the  original.  Again,  research 
activities  concentrated  on  the  pitch  extraction  algorithm.  The  advent 
of  serial  pitch  extraction  and  its  inclusion  and  success  in  the  trans¬ 
mitter  ordained  its  use  in  the  receiver. 

In  previously  investigated  algorithms,  the  pitch  extraction  algo¬ 
rithm  of  the  receiver  was  identical  to  that  of  the  transmitter.  An 
Important  innovation  of  this  final  algorithm  is  that  there  exists  a 
major  difference  in  the  receiver's  algorithm  as  compared  to  the  trans¬ 
mitter's.  Noting  that  the  receiver  has  a  reduced  amount  of  pitch  inf or 
mation  available  to  it,  the  receiver  algorithm  is  modified  to  more 
judiciously  use  this  information.  The  receiver  correlates  backward  as 
well  as  forward  one  searching  blocklength  for  pitch  information.  The 
transmitter  algorithm  only  looks  forward. 

The  product  of  the  receiver  is  a  sequence,  s(k),  which  is  the 
approximation  of  the  input  speech  sequence,  s(k).  The  rate  is  equal  to 
the  input  rate  of  8000  Hz.  The  similarity  s(k)  and  s(k)  is  used  as  an 
indication  of  the  performance  of  the  system.  Therefore,  side-by-side 
taped  comparisons  are  used  as  a  major  indication  of  performance. 


CHAPTER  IV 

INITIAL  ALGORITHM  REFINEMENTS 

4.1  Introduction 

The  initial  algorithm  refinements  are  described  in  this  chapter. 
These  refinements  are  those  that  are  principally  required  by  the  man¬ 
dates  of  Chapter  I.  The  hardware-induced  modifications  are  deferred  to 
Chapter  V.  Chapter  I  imposes  two  considerable  requirements.  It  re¬ 
quests  a  transmission  rate  of  16000  bits/second  and  a  bit  error  rate  of 
0.1Z.  These  conditions  immediately  translate  to  significant  algorithm 
adjustments.  Although  both  stipulations  are  channel  characteristics,  it 
is  necessary  to  modify  the  transmitter  and  receiver  algorithms  since  the 
channel  is  considered  unalterable. 

The  system,  as  it  appeared  following  these  initial  changes,  is 
illustrated  in  Fig.  4.1.1.  Hardware  induced  refinements  that  arose 
after  these  initial  changes,  such  as  the  switch  to  8  kHz  sampled  input 
speech,  are  considered  in  Chapter  V.  The  discussion  of  their  effects  on 
the  refinements  presented  in  this  chapter  are  also  deferred  to  Chapter 
V.  The  following  discussion  is  subdivided  into  two  sections.  The  first 
describes  the  refinements  due  to  the  new  transmission  rate.  The  second 
relates  the  added  alterations  caused  by  the  BER  condition. 

4.2  16  kb/s 

The  16  kb/s  transmission  rate  and  the  use  of  a  6400  Hz  sampled 
speech  signal  affected  the  ARC  and  inverse  ARC  considerably.  The  high 
transmission  rate  and  the  low  sampling  rate  allowed  the  use  of  a  fixed 
length  (non-run  length)  code.  A  bit  rate  of  16  kb/s  and  compressed 
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gure  4.1.1  Stage  of  Algorithm  Refinements 


speech  rate  of  3200  Hz  entitles  the  algorithm  to  an  average  of  5  bits 
per  sample.  A  five  bit  fixed  length  code  allows  for  a  32  level  quan¬ 
tizer.  Since  the  quantizer  is  symmetric  about  its  center,  an  odd  number 
of  quantizer  levels  are  used.  Therefore,  the  sampling  frequency  and 
transmission  rate  dictate  a  31  level  quantizer  if  a  fixed  code  is  used. 
Since  an  eleven  level  quantizer  produced  good  results  [Pande,  1981]  in 
the  original  algorithm  a  31  level  quantizer  is  initially  accepted  (and 
experimentally  supported  in  subsequent  tests). 

Earlier  Investigations  involved  fewer  levels  (11,  for  Instance). 

In  these  investigations,  the  thresholds,  scaling  factors,  and  expansion 
factors  are  judiciously  hand-picked.  (The  thresholds  and  scaling  fac¬ 
tors  are  the  input  and  output  values  of  the  quantizer  and  the  expansion 
factors  are  the  scalars  used  to  update  the  standard  deviation  computa¬ 
tion.)  The  advent  of  31  levels  brings  with  it  the  need  for  a  more 
systematic  method  of  choosing  these  values.  Because  the  magnitudes  of 
speech  samples  are  inherently  exponentially  distributed,  a  method  which 
Incorporated  a  logarithmic  scaling  of  the  quantizer  levels  was  pursued. 
The  final  system  allows  for  the  choice  of  a  base  and  a  maximum  scaling 
factor.  The  method  is  illustrated  in  Fig.  4.2.1.  In  this  figure,  the 
following  abbreviations  and  definitions  are  employed 

NQ  ■  number  of  quantization  levels, 

MSF  ■  maximum  scaling  factor, 

SFq  -  q*h  scaling  factor,  and 
THq  ■  q**1  threshold. 

The  quantizer  is  symmetric  about  zero  because  speech  is  nearly  as 
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Generated  Scaling  Factors  and  Thresholds 
DETERMINATION  OF  SCALING  FACTORS  AND  THRESHOLDS 


negative-going  as  it  is  positive-going.  Therefore,  the  zero-level  is 
pre-chosen  and  there  are  only  Q  Independent  levels,  where 

2 

The  scaling  factor  is  mapped  to  the  value  K  on  the  ordinate  by 

the  shifted  logarithmic  function 

Y  -  logioCX+1). 

(X  being  the  independent  variable.)  The  ordinate  is  then  equally 
divided  into  Q  sections  from  0  to  K  where 

R  -  logxoCMSF+1) 

The  ordinate  values  are,  therefore, 

0 

<1/Q)K 

<2/Q)R 

• 

e 

e 

«Q-2)/Q)K 

((Q-1)/Q)K 

K 

It  is  desired  to  map  these  ordinate  values  back  to  the  abscissa  to 
obtain  the  logarithmically  scaled  values.  To  add  flexibility  to  this 
technique,  the  base  of  the  logarithm  used  to  map  the  values  from  the 
ordinate  to  the  abscissa,  termed  the  distributing  base,  was  made 


The  effectiveness  of  this  technique  is  justified  by  a  comparison 
between  the  methodically  chosen  factors  and  the  hand-picked  factors.  A 
minimum  mean  square  error  is  used  to  compare  the  two  sets  of  data.  Other 
comparisons,  such  as  percent  error  were  considered  and  discarded.  Their 
undesirable  weighting  of  small  errors  in  the  lower  scaling  factors  led  to 
their  disqualification.  Using  the  mean  square  error  criterion,  a  local 
minuanim  mean  square  error  was  found  to  be  0.4  and  occurred  at  8  ■  9.0. 

As  in  proceeding  work  [Melsa  and  Pande,  1981],  the  thresholds  were 
chosen  to  be  equidistant  from  successive  scaling  factors. 

The  hand-picked  data  set  used  in  the  comparison  is  depicted  in  Fig. 
4.2.2  and  the  resulting  logarithmically  distributed  scaling  factors  are 
illustrated  in  Fig.  4.2.3.  Both  of  these  figures  incorporate  the  usual 
input-output  design  of  past  publications.  The  scaling  factors  and 
thresholds  agree  nearly  exactly  except  for  one  level  in  the  middle  of 
each  side.  This  abberatlon  is  momentarily  Ignored. 

The  expansion  factors  are  the  parameters  that  remain  unasslgned. 

The  expansion  factors  are  used  to  force  the  normalized  error  to  the 
middle  of  the  quantizer.  The  approximation  to  the  standard  deviation 
should,  therefore,  remain  unchanged  if  the  center  of  either  half  of  the 
quantizer  is  engaged.  For  this  reason  and  due  to  the  availability  of  a 
large  number  of  factors,  the  middle  three  factors  are  set  to  unity.  The 
remaining  upper  and  lower  values  are  distributed  linearly  between  the 
largest  and  smallast  expansion  factors,  respectively.  An  Altken  double 
sweep  multivariate  search  [Luenberger,  1965]  determined  the  local  maximum 
of  the  SNR;  pinpointing  the  maximizing  maximum  and  minimum  expansion 


values 


e  (k)/<r(k) 


0  PICKED  SCALING  FACTORS  AND  THRESHOLDS 


1  :j*iM**;st.#  *m*<£€} gg^  4  ,--•; 


A  table  listing  of  the  quantizer  thresholds,  scaling  factors,  and 
expansion  factors  for  speech  at  a  6400  Hz  and  a  8000  Hz  sampling  rate 
is  presented  in  Fig.  4.2.4.  The  large  number  of  expansion  factors 
allowed  a  tighter  tracking  of  the  normalized  error  sequence's  move¬ 
ments.  This  reduced  the  occupancy  frequency  in  the  outer  levels  of  the 
quantizer.  The  maximum  scaling  factor  was  appropriately  reduced  to 
improve  the  use  of  the  quantizer  levels. 

4.3  Bit  Error  Bate 

The  low  channel  bit  error  rate  of  0.1Z  precipitated  a  noiseless 
channel  approximation.  The  simulation  results,  therefore,  established 
a  tight  upper  bound  to  the  system  performance  in  the  presence  of  a  0.1% 
BER  channel.  The  assumption  of  no  channel  errors  facilitated  the  simu¬ 
lation  considerably . 

An  outcome  of  this  assumption  is  the  inclusion  of  the  coding  and 
decoding  blocks  in  the  channel  model.  Since  no  errors  are  generated, 
the  estimated  quantizer  levels  q(k),  of  the  receiver,  are  exact  replicas 
of  the  original  quantizer  levels,  q(k).  There  is  no  reason,  therefore, 
to  model  the  coding  and  decoding  blocks.  They  are  deemed  nonessential 
and  removed  from  the  transmitter  and  receiver.  These  blocks  are,  of 
course,  essential  to  a  hardware  system  and  for  this  reason,  they  are 
considered  as  part  of  the  channel. 

This  concept  is  illustrated  in  Fig.  4.3.1.  There,  the  bold  lines 
denote  the  simulation  flow  of  information.  The  dotted  lines  illustrate 
the  imagined  information  path  through  the  noiseless  channel. 
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Figure  4. 2. A  Quantizer  Thresholds ,  Scaling  Factors,  and  Expansion 
Factors  for  Speech  at  6400  Hz  and  8000  Hz  Sampling. 
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The  removal  of  the  buffers  In  the  transmitter  and  receiver  is 
assumed  when  this  modification  is  made.  The  buffers  are  removed  due  to 
a  combination  of  the  fixed  length  code  condition  and  the  temporary 
assumption  of  no  channel  errors.  The  removal  of  the  buffer  takes  with 
it  the  buffer  control  signal.  Since  there  is  no  error  protection  or 
correction  coding  needed  and  further  since  a  fixed  length  code  is  used, 
transmission  need  not  be  buffered.  The  original  system  incorporated  a 
variable  threshold  quantizer  that  adapted  according  to  a  buffer  control 
signal.  This  is  now  unnecessary. 

The  modifications  noted  in  this  chapter  are  justified  by  a  side-by- 
side  comparison  of  processed  and  unprocessed  speech  sentences  and  the 
predictor's  computed  performance.  At  this  point,  the  SNR  of  the  pre¬ 
dictor  was  recorded  in  excess  of  22  dB  for  sentence  1  and  over  27  dB  for 
sentence  6.  A  side-by-side  listening  session  established  the  speech 
estimate  as  being  a  very  good  representation  of  the  original  speech. 


Figure  4.3.1  Channel  with  Temporary  BER  Assumption 


CHAPTER  V 


HARDWARE  INDUCED  MODIFICATIONS 

5.1  Introduction 

Th«  refinements  that  are  principally  induced  by  hardware  considera¬ 
tions  are  presented  in  this  chapter.  It  is  subdivided  into  four  sec¬ 
tions.  They  are  1)  serial  pitch  picking,  2)  8  kHz  sampling  rate,  3) 

hardware  compatible  ARC,  and  4)  finite  word  length.  The  section  titles 
refer  to  the  hardware  desired  condition  which  instigated  the  investiga¬ 
tion  into  that  portion  of  the  algorithm. 

5.2  Serial  Pitch  Picking 

The  original  algorithm  chose  pitch  by  first  buffering  a  block  of 
speech  samples.  The  first  and  last  third  of  the  samples  in  the  block 
are  then  scanned  for  the  maximum  speech  sample  magnitude  and  the 
clipping  level  is  set  to  a  specified  percentage  of  this  maximum  value. 

L  2 

CK  -  y  *max{  S(n)  :  n  e  [  0,  —  ]  or  n  e  t  —  L,  L  ]} 

where  L  is  the  block  length  in  samples,  Cg  is  the  clipping  level  for  the 
k^h  determination  of  pitch,  and  s(n),  n  e  [0,L]  are  the  speech  samples 
in  the  buffer.  Another  buffer  of  the  same  length,  L,  is  filled  by  three 
value  center  clipping  the  buffered  speech  samples.  A  correlation  compu¬ 
tation  is  performed  on  the  clipped  speech  and  a  shifted  version  of  it. 
The  correlation  is  computed  over  the  specified  search  range  for  each 
shift  between  its  maximum  and  minimum  shifts.  The  shift  value  associat¬ 
ed  with  the  maximum  correlation  computation  is  accepted  as  the  pitch 
period.  A  window  is  formed  using  this  value  and  compression  or 


expansion  Is  performed.  This  technique  of  pitch  detection  is  especially 
attractive  due  to  its  simple  hardware  implementation. 

After  the  transmitter  compresses  the  speech,  the  speech  buffer  is 
shifted  by  two  pitch  periods  and  refilled  with  new  samples.  By  shifting 
out  two  pitch  periods  of  processed  speech  and  shifting  in  two  pitch  per¬ 
iods  of  new  speech,  the  next  clipping  level,  Cg+1,  can  be  different  than 
the  previous  clipping  level,  Cg.  This  is  caused  by  the  term  "max{»}"  in 
the  clipping  level  computation,  which  can  result  in  different  values  for 
the  two  speech  sample  buffers.  Since  the  entire  buffer  of  speech  is 
clipped  using  a  single  clipping  level,  the  unprocessed  samples  of  the 
previous  speech  buffer  must  be  reclipped  with  the  new  clipping  level. 
This  means  that  the  entire  previously  clipped  speech  is  useless  and  must 
be  discarded.  Furthermore,  clipping  cannot  take  place  until  the  speech 
buffer  is  completely  filled  and  a  new  dipping  level  is  determined. 

The  problem  of  waiting  for  the  speech  buffer  to  fill  before  start¬ 
ing  the  clipping  process  and  the  problem  of  clipping  speech  more  than 
once  before  processing  the  associated  speech  samples  are  eliminated  by 
serial  pitch  picking  [Arora,  1983].  The  solution  to  this  time  wasting 
problem  is  termed  serial  pitch  picking.  The  serial  technique  uses  a 
scaled  approximation  to  the  speech  signal  energy  as  a  sample  by  sample 
updated  dipping  level.  The  clipping  level,  Cg,  is 

Cfc  ■  (J»Eg 

where 

Eg  -  a»Eg_i  +  (l-rs)|Sg|. 

8  is  the  scaling  constant  and  a  is  a  constant  related  to  the  length  of 


the  transient  response.  The  response  of  this  computation  is  described 

by 

H(Z)  -  ( l-o )  ~- 
Z-u 

which  describes  a  low  pass  filter  whose  time  constant,  in  samples,  is  a 
function  of  a.  The  time  constant,  n,  in  samples  is 

~  1 

n  * -  . 

1-a 

As  the  new  speech  enters  the  buffer,  the  energy  is  updated,  the 
clipping  level  is  determined,  and  the  speech  is  three  value  center 
clipped  using  this  value.  The  clipped  speech  sample  is  shifted  into  the 
clipped  speech  buffer  and  the  routine  is  ready  to  accept  the  next  speech 
sample.  If  the  previously  detected  pitch  period  is  large,  a  large  num¬ 
ber  of  samples  must  be  accepted  and  clipped  to  set  up  the  buffers  for 
the  next  pitch  detection.  The  large  number  of  new  samples  will  take  a 
proportionately  large  hardware  time  period  to  clip.  This  will  not  pres¬ 
ent  a  problem  since  the  time  between  pitch  detections  is  also  propor¬ 
tional  to  the  pitch  period  length.  (As  the  number  of  samples  Involved 
in  the  compression  increases,  the  time  required  to  compress  these  sam¬ 
ples  increases  and  the  frequency  of  pitch  detection  decreases.) 

The  unprocessed  speech  of  the  previous  speech  buffer  was  clipped  at 
the  correct  clipping  levels.  The  clipped  speech  associated  with  the  un¬ 
processed  speech  is,  therefore,  retained  and  the  redundancy  of  the  clip¬ 
ping  in  the  block  pitch  picking  is  eliminated.  The  clipped  speech  is 
generated  as  the  samples  enter  the  speech  buffer.  Therefore,  the 


problem  of  waiting  for  the  speech  buffer  to  fill  found  in  block  pitch 
picking  is  also  eliminated. 

The  computational  simplicity  of  the  three  value  center  clipped 
speech  correlator  is  sustained.  Therefore,  serial  pitch  picking  retains 
the  simple  hardware  implementation  of  the  block  pitch  picking  and  yet 
eliminates  the  problems  associated  with  it.  The  comparison  of  the  two 
techniques  is  illustrated  in  Fig.  5.2.1.  The  difference  in  the  algorithm 
for  the  receiver  is  that  at  most  one  maximum  pitch  period  of  samples 
could  be  used  whereas  the  transmitter's  maximum  is  twice  this  value. 

Numerous  tests  were  performed  and  the  suggested  parameter  values 

are 

0  -  1.5 
o  -  0.992 

This  a  value  allowed  the  proper  tracking  of  the  signal  energy  required  by 
the  technique  and  the  above  0  value  insures  the  success  of  the  clipping 
routine  by  properly  scaling  the  resultant  signal  energy.  The  fundamental 
concept  of  three  value  center  clipping  is  to  retain  only  the  first  couple 
of  positive  and  negative  going  pulses  of  each  pitch  period  [Schafer, 
1976].  Figure  5.2.2(a)  shows  how  block  pitch  extraction  succeeds  in 
clipping  in  this  manner.  Figure  5.2.2(b)  demonstrates  serial  clipping's 
similar  performance  and  Fig.  5.2.2(c)  is  the  independent  result  of  the 
two  methods. 

Two  problems  developed  in  the  serial  clipping  implementation.  Half 
pitch  was  sporadically  picked  in  the  voiced  segments  and  a  buzzing  sound 
appeared  in  the  output  speech  during  unvoiced  segments.  This  buzzing 
sound  is  attributed  to  the  repeated  minimum  pitch  periods  picked  in 
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unvoiced  segments  by  Che  seriel  pitch  extrectlon  technique.  Both  prob¬ 
lems  ere  cherged  to  the  clipped  epeech  sequence,  which  alone ,  affects 
the  correlation  coaputatlon. 

The  half-pitch  problem  is  solved  by  implementing  a  digital  low  pass 
filter  and  by  modifying  the  processing  of  the  correlation  computations. 
The  actual  correlation  computation.  Itself,  remained  unaltered.  How¬ 
ever,  the  processing  of  the  elements  of  the  correlation  vector  was 
changed.  The  routine  now  seeks  the  largest  of  the  correlation  vector's 
one  sample  local  maxima.  A  one  sample  local  maximum  is  defined  as  an 
element  whose  magnitude  is  greater  then  or  equal  to  its  immediate  neigh¬ 
bors.  The  largest  local  maxima  indicates  the  pitch  period  as  did  the 
vector's  maximum  element  in  the  original  algorithm.  The  incorporation 
of  this  strategy  eliminated  the  possibility  of  choosing  either  the  maxi¬ 
mum  or  minimum  pitch  when  they  represent  the  rising  or  falling  side, 
respectively,  of  a  correlation  sequence's  local  maximum.  This  technique 
aided  in  picking  true  pitch  in  the  voiced  regions. 

The  half-pitch  problem  is  further  mastered  by  installing  a  low  pass 
filter  between  the  speech  sequence  and  the  clipper.  It  is  a  9th  order 
linear  phase  FIR  filter  (see  Appendix  D).  The  filter  is  designed  to 
smooth  out  the  anomalies  of  the  individual  voiced  segments,  thereby 
forcing  the  pitch  periods  to  appear  more  similar.  This  aids  the  corre¬ 
lator  by  coercing  the  clipped  speech  of  each  pitch  period  to  be  more 
alike.  The  result  of  these  two  methods  is  the  elimination  of  half -pitch 
picked  in  voiced  segments. 

The  problem  of  repeated  minimum  pitch  period  picked  in  the  unvoiced 
segments  is  corrected  by  implementing  a  controlled  random  pitch 


generator.  Whenever  the  minimum  pitch  period  is  picked,  a  pseudo-random 
pitch  is  chosen  and  employed  as  the  correct  pitch.  The  random  pitch  is 
restricted  to  the  lower  half  of  the  search  range  so  that  a  large  random¬ 
ly  picked  pitch  period  in  an  unvoiced  region  will  not  adversely  affect 
the  beginning  of  a  voiced  segment.  Since  the  pitch  in  the  voiced 
regions  is  considerably  larger  than  the  minimum  value,  the  random  pitch 
picking  should  not  be  engaged  and  no  effect  to  the  voiced  regions  should 
occur.  If  a  minimum  pitch  period  is  picked  in  a  voiced  region,  it  is  an 
incorrect  pitch  and  randomly  choosing  another  pitch  period  will  not  harm 
and  may  improve  the  system  performance.  This  variation  to  the  pitch 
picking  algorithm  eliminates  repeated  minimum  pitch  estimations  in  the 
unvoiced  regions.  The  problem  is,  therefore,  solved. 

For  the  first  time,  the  receiver  pitch  extraction  algorithm  differed 
from  that  of  the  transmitter.  It  is  noted  that  a  forward  searching  cor¬ 
relator  in  the  receiver  has  less  pitch  information  available  to  it  than 
did  its  counterpart  in  the  transmitter.  The  larger  probability  of  in¬ 
correctly  choosing  pitch  at  the  boundaries  of  voiced  regions  results  in 
serious  distortions  of  the  output  speech.  Therefore,  whenever  pitch  is 
present  in  the  compressed  signal  estimate  y(k),  it  is  critical  that  the 
information  is  extracted  and  used.  The  correlator  has  been  modified  to 
Include  a  larger  searching  range.  The  beginning  of  this  search  range  is 
located  in  previously  expanded  compressed  speech.  The  buffer  sizes  in 
the  receiver  must,  therefore,  be  larger  since  this  previously  processed 
speech  is  normally  discarded.  This  variation  led  to  Improved  performance 
of  the  previously  unsatisfactory  receiver  serial  pitch  extraction  algo¬ 
rithm.  Further  testing  into  the  possibilities  of  reducing  the  search 


range  from  20  to  100  to  a  digital  hardware  compatible  range  Sueh  as  24 
to  88  was  performed.  Multiple  testing  on  the  boundaries  of  the  search¬ 
ing  range  of  the  correlator  proved  that  no  perceptable  difference  in 
speech  quality  accompanied  the  reduction  of  the  searching  range  as  far 
as  24  to  88.  This  searching  range  is.  therefore,  adopted  in  the  final 
algorithm. 

5.3  8  kHz  Sampling  Rate 

Throughout  the  term  of  this  research,  studies  paralleling  those 
described  for  a  system  operating  on  a  6400  Hz  sampled  speech  signal  were 
performed  on  an  alternate  system  designated  for  a  8000  Hz  sampling  rate. 
The  8  kHz  sampling  rate  is  a  standard  industry  rate  and  repeatedly  pre¬ 
sented  itself  as  a  possible  alternative  to  the  6400  Hz  system.  Eventu¬ 
ally,  the  superior  availability  of  reliable  8  kHz  filters  proved  to  be 
the  convincing  argument  for  the  shift  to  an  8  kHz  system. 

The  modifications  necessary  Included  a  reduction  in  the  number  of 
quantization  levels  to  -1  ■  15,  instead  of  the  31  levels  previously 
used.  The  31  levels  of  the  6.4  kHz  system  produced  such  fine  quantiza¬ 
tion  that  a  fixed  (first  order)  linear  predictor  was  considered.  The 
quality  of  the  fixed  predictor  system  at  6.4  kHz  is  identical  to  that  of 
the  adaptive  predictor  at  6.4  kHz.  The  reduction  in  available  quantiza¬ 
tion  levels  and  the  goal  to  produce  equivalent  quality  speech  at  an  8 
kHz  sampling  rate,  dictated  the  use  of  an  adaptive  predictor.  The 
thresholds,  scaling  factors,  and  expansion  factors  are  all  selected  us¬ 
ing  the  technique  that  generated  the  31  level  values  presented  in 
Chapter  IV. 


Other  parameters  that  require  adjustment  are  the  alpha  used  in  the 
signal  energy,  the  block  length  and  searching  range  of  the  correlator, 
and  the  LPF  coefficients.  The  alpha  value  reported  in  the  previous  sec¬ 
tion  is  the  corrected  value  for  the  8  kHz  system.  The  searching  range 
limits  are  changed  to  40  to  110,  which  reflects  nearly  an  8000/6400  * 

5/4  scaling  factor.  The  8  kHz  filter  with  1000  Hz  cutoff  freqency  pre¬ 
sented  in  the  Appendix  is  the  filter  used  in  the  8  kHz  sampling  rate 
system. 

5.4  Hardware  Compatible  ARC 

The  main  processing  Integrated  circuits  of  the  proposed  hardware 
implementation  are  the  INTEL  8085  and  NEC  7720.  Although  each  in  itself 
is  a  very  versatile  chip,  neither  supports  division  capabilities.  This 
hardware  consideration  leads  to  significant  modifications  in  the  adap¬ 
tive  residual  coder.  Specifically,  the  adaptation  of  predictor  coeffi¬ 
cients  and  the  normalization  of  the  input  to  the  quantizer  use  a  divi¬ 
sion  operation  which  must  be  eliminated. 

The  original  system  adapted  the  coefficients  according  to  the 
steepest  descent  of  e2(k)  [Melsa  and  Cohn,  1975]. 

atCk+l)  -  5bi  +  (ln5)[ai(k)  + 

<|y(k)|>2 

where  <|y(k)|>  is  a  biassed  exponential  time  average  of  |y(k)|  [Melsa 
and  Panda,  1981]. 

The  division  by  <|y(k)|>2  must  be  modified.  This  is  accomplished 


by  adopting  a  new  adaptation  procedure.  The  final  algorithm  uses  a 


fixed  gain  (percentage)  adaptation. 

*i(k)  -  b4  +  ai*(k) 

»i*(k+l)  -  3(aj*(k)  +  sgn{e(k),y(k-i)}  max{g|ai*(k)  |  ,at,iag}) 
at(k+l)  -  bi  +  at*(k+l) 

where  the  recommended  values  are 

bi  -  1.79 

b2  -  -1.40 
b3  -  0.57 

b4  -  -0.15 
3  -  0.95 

ablas  "  0.001 
g  -  0.15 

The  b*  are  determined  using  the  technique  suggested  by  R.A.  McDonald 
(1966).  The  term,  max{«} ,  computes  the  biased  percent  change  desired. 
The  a^ -tas  term  helps  the  coefficient  move  from  the  zero  level.  Without 
such  a  term,  a  coefficient  which  attained  a  zero  value  would  remain  at 
zero.  The  term,  sgn{*} ,  is  a  first  order  estimate  of  the  direction  in 
which  the  coefficient  should  move  to  improve  the  predictor  performance. 

The  original  quantizer  accepted  the  error  sequence  from  the  predic¬ 
tor  and  normalized  these  values.  The  normalized  values  were  then  com¬ 
pared  to  the  thresholds  and  mapped  to  appropriate  quantization  levels, 

•(k)H0FM  "  «(k)/o(k) 

•OOHOSM  a  T(i)  -  SF(i) 


The  division  by  the  approximation  to  the  standard  deviation,  <j(k),  must 
be  eliminated.  This  division  is  removed  by  comparing  the  error 


sequence,  e(1c),  to  adaptive  thresholds, 

e(k)  a  T(i,k)  -  SF(i,k) 

T(i,k)  -  o(k)  x  T(i) 

SF(i,k)  -  o(k)  x  SF(i) 

These  modifications  resulted  in  no  noticeable  difference  in  output 
speech  quality. 

5.5  Finite  Word  Length 

With  any  real  time  system  a  restriction  on  the  word  length  is 
imposed.  The  integrated  circuitry  proposed  for  the  construction  of  a 
real  time  system  uses  16  bit  words.  Therefore,  a  software  simulation 
which  restricts  all  wordlengths  to  a  maximum  of  16  bits  is  desired.  It 
should  be  noted  that  it  is  possible  to  Increase  the  word  accuracy  to 
some  integer  multiple  of  16,  such  as  32.  This  requires  a  much  more 
complicated  hardware  algorithm  since  each  software  computation  would  in 
fact  require  more  than  one  hardware  computation.  Keeping  track  of  the 
order  and  location  of  multiple  word  computations  further  complicates 
such  a  scheme.  Therefore,  a  16-bit  finite  word  length  simulation  is 
highly  desired. 

In  fact,  a  16-bit  finite  wordlength  simulation  resulted  in  very 
good  quality  output  speech  as  compared  to  the  output  of  the  simulation 
having  no  such  restriction.  In  all  computations  and  representations  the 
most  significant  bit  is  left  unoccupied  to  give  a  factor  of  2  margin  of 


error.  This  is  done  to  virtually  eliminate  the  possibility  of  severely 
distorting  overflows.  Reiterating,  the  performance  of  the  finite  word- 
length  simulation  is  comparable  to  that  of  the  final  unrestricted  word 


length  simulation.  A  16-bit  wordlength  simulation  is,  for  this  reason, 
adopted. 

The  speech  quality  as  perceived  by  side-by-side  comparisons  of  a 
taped  output  version  of  the  speech  and  input  was  qualified  as  very  good 
The  predictor  SNR  is  in  excess  of  23  dB  for  sentence  1. 

1 


40 


v 


JtNJ‘ 


CHAPTER  VI 


METHODOLOGY 

6.1  Introduction 

This  chapter  introduces  the  software  capabilities  and  describes  the 
associated  output  files  of  the  simulation  in  the  context  of  their  use¬ 
fulness  through  past  research  activities  and  for  the  benefit  of  future 
researchers  using  this  software.  The  output  files  are  presented  first 
to  simplify  the  task  of  explaining  the  program  options  which  employ 
these  files.  They  are  Illustrated  as  output  files  in  Fig.  6.1.1*  The 
file  muses  will  precede  their  Individual  discussion  and  a  (T)(R)  at  the 
beginning  of  a  file  name  indicates  that  an  option  exists  to  output  this 
file  from  the  transmitter  and/or  the  receiver.  If  it  Is  generated  in 
the  transmitter,  the  first  letter  of  the  name  is  a  T,  and  if  the  receiv¬ 
er  is  the  source ,  the  file  name  begins  with  an  R. 

6.2  Files  and  Options 

The  changes  Incorporated  in  the  pitch  extraction  algorithms  result¬ 
ed  in  the  generation  of  a  number  of  files  which  aided  in  the  evaluation 
of  the  technique  under  consideration.  Those  associated  with  this  topic 
are  (T)(R)  PITCH.dat,  (T)(R)  CLPSP.dat,  (T)(R)  SENGY.dat  and  (T)(R) 
FILTER.dat.  (T)(R)  PITCH.dat  proved  to  be  a  vital  file  in  the  investi¬ 
gation  of  pitch  extraction.  It  contains  the  values  of  pitch  picked  in 
the  respective  algorithm,  transmitter  or  receiver.  A  simple  listing  of 
these  files  contains  a  large  quantity  of  Information.  With  this  list¬ 
ing,  it  is  possible  to  make  a  quick  determination  of  whether  the  pitch 
is  being  extracted  properly  or  not.  Given  a  sampling  rate  of  8000  Hz, 
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the  expected  pitch  period  length  in  voiced  male  segments  is  approximate¬ 
ly  60  samples.  The  algorithm’s  agreement  with  this  can  be  expeditiously 
evaluated  from  the  listing.  A  critical  condition  of  the  pitch  extrac¬ 
tion  algorithm  is  that  it  should  choose  the  same  pitch  period  in  the  re¬ 
ceiver  as  it  chose  in  the  transmitter.  Side-by-side  comparison  of  the 
files  will  determine  if  this  condition  is  fulfilled. 

PITCH.dat  illustrated  the  periodicity  of  pitch  picked  in  silence 
regions  by  the  serial  pitch  extraction  algorithm  and  precipitated  the 
implementation  of  the  random  pitch  picker.  PITCH.dat  forms  the  input 
of  the  program  titled  PICHMAK.ftn  which  transforms  the  pitch  periods 
into  a  series  of  rectangular  pulses  of  height  equal  to  the  pitch  period 
and  duration  equal  to  twice  the  pitch  period.  The  output  of  this  file 
is  compatible  input  for  the  plotting  routine  used  throughout  the  re¬ 
search  called  DATAPLOT.ftn.  The  usefulness  of  the  hard  copy  output  is 
that  it  can  be  laid  upon  a  speech  segment  giving  a  pictorial  account  of 
the  pitch  extraction  algorithm's  performance.  This  is  Illustrated  in 
Fig.  6.2.1.  Although  this  is  a  more  precise  and  reliable  report  of  the 
effectiveness  of  the  pitch  extraction  routine,  the  cost,  in  time,  re¬ 
quired  to  8 pool  the  results  to  the  plotter,  necessitates  the  use  of  the 
pitch  period  lists  in  multiple  run  testing. 

(T)(R)CLPSP.dat  is  a  list  of  the  three  value  center  clipped  speech. 
It  is  also  compatible  with  the  plotting  program,  the  output  of  which  de¬ 
termined  the  potency  of  the  serial  clipping  routine.  The  expected  out¬ 
put  of  such  a  clipping  routine,  with  speech  as  an  input,  is  a  series  of 
positive  one's  occurring  at  the  large  positive  going  section  of  the 
pitch  pulses  and  negative  one's  for  their  negative  going  counter  parts. 
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Zeros  are  expected  to  fill  in  the  gaps.  Laying  this  on  top  of  the  input 
speech  file  is,  again,  effective  in  ascertaining  the  clipping  action. 

The  tine  associated  with  obtaining  the  hard  copy  should  be  evaluated  be¬ 
fore  performing  a  multiple  run  evaluation.  This  file  illustrated  the 
proximal  relationship  between  the  block  pitch  extraction's  clipping 
routine  and  that  of  the  serial  pitch  extractor. 

(T)(R) SENGY.dat  is  a  listing  of  filtered  speech  signal  energy  (com¬ 
pressed  speech,  in  the  receiver).  Used  with  the  plotting  program,  this 
file  demonstrated  the  effectiveness  of  the  signal  energy  computation. 

It  is  necessary  to  compare  the  signal  energy  plot  to  the  output  of  the 
filter  and  not  to  the  original  speech,  since  clipping  will  be  performed 
on  the  filtered  speech.  Comparison  of  the  signal  energy  to  the  filtered 
speech  assisted  in  calibrating  the  clipping  value,  CLPP.  CLPP  is  the 
scaling  value  of  the  signal  energy;  the  product  being  the  instantaneous 
clipping  level. 

(T) (R)FILTER.dat  is  that  file  whose  hard  copy  is  compared  to  the 
signal  energy  computation  as  discussed  in  the  previous  paragraph.  It  is 
the  low  pass  filter  output  of  the  pitch  extraction  units.  It's  plot 
verified  the  filter's  smoothing  effect  on  the  input  speech. 

The  combination  of  these  files  depicted  the  performance  of  each 
component  of  the  pitch  extraction  unit.  Their  aid  was  unparalleled  in 
the  modification  of  the  TDHS  algorithms  and  it  is  suggested  that  they 
continue  to  be  used  in  the  evaluation  of  any  further  simulation  adjust¬ 
ments. 

TZHAT.dat.  (simply,  ZHAT.dat .  in  the  transmitter)  contain  the  com¬ 
pressed  speech  files.  A  side-by-side  listening  session  comparing  a 
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taped  version  of  these  files,  as  well  as,  the  examination  of  their  hard 
copy  plots  greatly  contributed  in  determining  the  potency  of  the  TDHS 
units.  The  hard  copy  is  expected  to  be  a  plot  highly  related  to  the 
original,  whose  time  axis  is  scaled  according  to  the  compression  ratio 
adopted.  The  taped  version  should  reflect  this  same  time  axis  scaling. 

QUANTjdat  and  DECO.dat  are  the  input  and  output  of  the  channel,  re¬ 
spectively.  They  are  sequences  of  quantization  levels  and  are  equivalent 

i 

when  the  noiseless  channel  approximation  is  engaged.  When  a  noisy  chan¬ 
nel  is  employed,  DECO.dat  is  the  error  infected  version  of  QUANT.dat. 

TSTATS.dat  is  a  file  output  by  the' transmitter  only.  It  contains  a 
presentation  of  the  system  performance  criteria  employed  in  the  original 
system.  TSTATS.dat  includes  items  such  as  the  average  signal-to-noise 
ratio  of  the  predictor  and  the  parameters  used  in  the  quantizer.  A 
sample  of  this  file  is  contained  in  Appendix  E.  Since  the  different  SNR 
measurements  pertain  to  the  transmitter  unit,  this  file  was  primarily 
commissioned  in  the  ARC  evaluations. 

SHAT.dat  contains  the  simulation's  estimate  of  the  original  speech 
file.  The  taped  reproduction  of  this  formed  the  most  highly  weighted 
estimation  of  the  system's  effectiveness.  Frequently,  the  quality  of 
the  output  file,  as  perceived  by  listeners,  was  the  only  criteria  used 
in  decisions.  This  file,  therefore,  formed  the  most  important  measure  of 
the  system's  performance. 

The  program  options  are  illustrated  in  Fig.  6.2.2.  There  are  four 
major  options  available  to  the  user.  The  simulations  available  are  an 
algorithm  that  exercises  the  TDHS  unit  only  (without  the  ARC)  and  one 
that  concatenates  the  TDHS  and  ARC  blocks.  Each  of  these  options  has  the 
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option  of  passing  the  pitch  picked  In  the  transmitter,  as  side  Informa¬ 
tion,  to  the  receiver.  No  cost  to  channel  capacity  was  evaluated  for 
the  separate  options  with  the  exception  of  the  option  of  TDHS-ARC  with 
no  pitch  passed.  This  option  represents  the  final  algorithm.  The  sepa¬ 
rate  options  assisted  In  assessing  the  noise  due  to  the  TDHS  algorithm, 
the  ARC  blocks,  and  the  Incorrect  pitch  in  the  receiver.  The  compari¬ 
sons  of  the  different  options  demonstrated  the  magnitude  of  each  of 
these  noise  sources. 

The  source  files  for  the  unrestricted  word  length  simulations  are 
presented  In  Appendix  A.  Reference  should  be  made  to  these  when  deciph¬ 
ering  Appendix  B,  the  finite  word  length  simulations ’  source  files. 

Both  sets  of  source  files  are  driven  and  fused  by  the  Included  Indirect 
command  file.  Appendix  C  contains  the  source  files  of  supporting  pro¬ 
grams  which  aided  In  the  development  process. 


CHAPTER  VII 


CONCLUSIONS 

The  final  system  Incorporates  all  of  the  specifications  of  the 
Chapter  I  and  fora*  a  valid  simulation  of  a  finite  word  length  time 
domain  harmonic  scaling-adaptive  residual  coder  system.  The  final  pre¬ 
dictor  SNR  is  in  excess  of  22  dB  for  the  sentence 

"Cats  and  dogs  each  hate  the  other." 

The  critical  topics  for  those  who  are  familiar  with  TDHS-ARC  system 
are  contained  in  Chapters  III,  IV,  and  V. 

The  issue  herein  reported  still  maintains  questions  unanswered  and 
forms  possible  topics  of  future  investigation.  The  optimization  of  the 
technique  used  to  distribute  scaling  factors ,  thresholds ,  and  expansion 
factors  is  one  example.  It  seems  that  the  maxima  of  these  sets  of 
parameters  could  be  reduced  to  the  extent  of  precipitating  added  accura¬ 
cy  in  the  finite  word  length  simulations.  An  investigation  into  the  re¬ 
moval  of  the  random  pitch  generator  is  pressed  by  the  effectiveness  of 
the  modified  correlator.  Random  pitch  is  picked  few  times  in  the  trans¬ 
mitter,  usually  at  the  beginning  or  end  of  the  sentence  and  rarely  is 
random  pitch  picked  successively.  The  receiver  results  are  similar. 

For  the  above  sentence,  the  receiver  picked  successive  random  pitch  only 
once.  It  does  not  seem  to  be  worthy  of  the  added  software  complexity  to 
implement  this  technique.  Investigations  into  the  difference  of  param¬ 
eter  values  of  the  receiver  as  compared  to  the  transmitter  need  to  be 
accomplished  (such  as,  the  time  constant  of  the  clipping  signal  energy 
computation) . 


The  final  algorithm  incorporates  the  suggested  TDHS-ARC  structure 
with  modifications  to  the  original  generated  by  the  elements  of  the 
original  proposal  and  by  hardware  restrictions.  For  the  TDHS  modifica¬ 
tions,  emphasis  was  placed  on  finding  a  hardware  compatible  pitch 
detector.  The  new  adaptation  scheme  is  the  major  variation  found  in  the 
ARC.  The  final  hardware  simulations  resulted  in  good  quality  speech  in 
the  presence  of  up  to  0.1Z  BER,  ascertained  by  side-by-side  comparison 
of  a  taped  reproduction. 

The  goal,  to  design  a  simulation  by  means  of  algorithm  refinements 
that  is  easy  to  implement  in  hardware  and  yet  does  not  degrade  the 
present  performance,  has  been  met. 


APPENDIX  A 


FINITE  WORD  LENGTH  SIMULATION  SOURCE  FILES 
The  necessary  nodules  to  run  a  simulation  are: 

1)  CMD16.CMD  (Indirect  command  file) 

2)  PFWL.FTN  (Parameter  builder) 

3)  TFWL.FTN  (Transmitter) 

4)  C8K.FTN  (Channel) 

5)  RFWL.FTN  (Receiver) 

The  task  transmitter  and  receiver  modules  should  be  built  from  the 
source  files  using  the  following  switches: 

>  UNITS  -  15 

>  ACTFILE  -  15 

All  other  source  files  use  the  default  values.  A  specific  file  manage 
ment  procedure  Is  dictated  by  the  indirect  command  file.  This  can 
easily  be  modified  to  fit  the  specific  needs  of  a  multiple  run  use. 
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TYPE  WHICH  PROGRAM  PATH  DO  YOU  WANT  TO  TAKE7' 

TYPE  *,*  DTDHS-ARC  WITH  PITCH  PASSED  TO  THE  RCVR* 

TYPE  *,*  2 )TOHS-ARC  WITH  PITCH  INDEPENDENTLY  GENERATED* 

TYPE  *,*  IN  THE  RCVR* 

TvpE  *.*  3  )TDHS  ONLY  WITH  PITCH  PASSED  TO  THE  TDHE  ALGORITHM' 
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APPENDIX  D 


FILTER  DESIGN 

FINITE  IMPULSE  RESPONSE  (FIR) , 
LINEAR  PHASE,  9TH  ORDER, 

DIGITAL  LOW  PASS  FILTER  SPECIFICATIONS 
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Wc  -  300  Hz 
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Wc  -  500  Hz 

h(l) 
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0.056 
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•  h(8)  ■ 
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0.021 
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-  h(7)  - 
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0.234 

0.025 

h<4) 

■  h(6)  » 

0.122 
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0.319 

0.028 

h(5) 

-  M5)  - 

0.130 

0.015 

0.336 

0.029 

8k//lk 

8k//. 5k 

6.4//lk 

6.4//. 5k 

Band  1: 

Lower  Band  Edge 

0.0 

0.0 

0.0 

0.0 

Upper  Band  Edge 

0.125 

0.063 

0.156 

0.078 

Band  2: 

Lower  Band  Edge 

0.131 

0.069 

0.164 

0.086 

Upper  Band  Edge 

0.500 

0.500 

0.500 

0.500 

APPENDIX  E 
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